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Abstract. The Griffiths inequalities for Ising spin-glass models with Gaussian 
randomness of non-vanishing mean are proved using properties of the Gaussian 
distribution and gauge symmetry of the system. These inequalities imply that 
correlation functions are non-negative and monotonic along the Nishimori line in the 
phase diagram. From this result, the existence of thermodynamic limit for correlation 
functions and pressure is proved under free and fixed boundary conditions. Relations 
between the location of multicritical points are also derived for different lattices. 



1. Introduction 

Two Griffiths inequalities provide a significant insight about the phase transitions in 
ferromagnetic Ising models pp. One of their formulations states that if the sets of 
Ising spins Sa = HjeA'S'i are coupled by an energy —JaSa (for some positive Ja), the 
free energy F and all the correlations are monotonic functions of the strength of any 
interactions J's, namely: 



These inequalities can be used in ferromagnetic Ising models to prove that the free 
energy and correlation functions have the thermodynamic limit under several boundary 
conditions and to demonstrate the existence of phase transitions for various lattices. 
The proof of equations (0) and Q [21 assumes that all the interactions among spins 
are ferromagnetic, a condition that clearly fails in spin glass models which have both 
ferromagnetic and antiferromagnetic interactions. 

Attempts to extend even only partially those or similar inequalities to the spin 
glasses have been unsuccessful until very recently after the technique of integration by 
parts has been powerfully exploited in the mathematically rigorous approaches: in 13] 
the main results are based on correlation inequalities that represents the equivalent of 
the first Griffiths inequality (0) for the Sherrington-Kirkpatrick model and the Edwards- 
Anderson model respectively. The monotonicity properties in [3 11] are proved not with 




(1) 



(2) 
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respect to the strength of the interaction but with respect to the variance of the random 
interaction. 

In this paper we show that a Gaussian Ising spin-glass model does fulfil the first and 
the second Griffiths inequalities (H)) and © with respect to the mean of the distribution: 
pressure and correlations are monotonic functions with respect to the mean. The results 
are proved on the Nishimori line (NL), a restricted space of the phase diagram in which 
several exact results follow from the gauge symmetries of the system jHllZI. 

We apply the resulting inequalities to prove the existence of the thermodynamic 
limit for the pressure and correlation functions under free and fixed boundary conditions. 
Moreover, we derive inequalities on the location of multicritical points. 

In the next section, we present our results and their proofs. The applications of 
the inequalities are discussed in the last section. Some of the details of calculations are 
described in the appendix. 



2. Inequalities 

Let us consider a finite box Q, subset of a regular lattice C. To each point i & Q we 
associate an Ising spin {Si = ±1) and denote their products by 

= n (3) 

We consider the spin glass model defined by the random potential 

Un=Y. PaJaSa , (4) 
Acn 

where /3a > is the inverse of local temperature of subset A and Ja is a quenched 
random variable which follows the Gaussian distribution with positive mean [Ja] = Jao 
and variance [{J a — ^ao)^] = The pressure function is defined as 

P=[log5:e^"], (5) 

s 

from which all physical quantities can be derived. 

The NL is defined in terms of the parameters > (for all the A) as [HI IZj 

(3a = — , Jao = o-axa ■ (6) 

O-A 

Our results are the following inequalities: 
dP 



dXB 



xb[{Sb + 1)]>0 (7) 



-[{Sc)] = 2xb[{{SbSc) - {SB){Scm > . (8) 



dxB 

Both inequalities hold for arbitrary subsets B,C as long as the parameters satisfy NL 
condition (jH)). 
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The first and second inequalities are proved using some properties of the Gaussian 
distribution and the gauge theory. To prove the first inequahty ((Tj), we observe that on 
the NL the formula of the total derivative gives 



dP dP dpB ^ dP dJeo 1 dP 



dP 



dxB OPb dxB ' OJbo dxB ctb d^B ' ' 
As shown in the appendix, we can derive the following identities from the properties of 
the Gaussian distribution for any operator O: 

dO 



[JBO] = JBo[0] + a 



B 



dJ, 



B 



d 



dJ, 



BO 



-[O] 



dO 



(10) 



Note that both equations (fTUj) and are valid on and away from the NL. The identity 
(HH) yields 

dP 



OJbo 

and a direct computation provides 
dP 



dp. 



B 



\Jb{Sb)] 



(12) 



(13) 



Using the gauge theory we have on the NL jHl EI 

[{JbSb)] = Jbo = ctbXb ■ (14) 

From equations (jH)), |T2|l - |T4|l . we immediately find equation ((Zj). 

It is also possible to obtain a related inequality, an explicit bound for the correlation 
function, as 



KSb)]>j^ {Be {A}) 



[{Sb)]>0 iB^{A}). 



For this purpose we observe that, when B G {^4}, 
[{JbSb)] -- 



n 

Acn 



Pa{Ja)Jb] [^Pa{Ja){Sb) 



(15) 



(16) 



Let us square both sides of the above equation and apply the Cauchy-Schwarz inequality 
to obtain 

[{JBSB)f<{crl + jl,)[{SB?]. (17) 

Since the gauge theory yields an identity ^ [7j 

[{Sb)] = [{Sb?] (18) 

on the NL under certain boundary conditions (free, periodic or fixed (all spins up)), we 
obtain equation (fT3jl from equations (fT^ and (fT7|) . When B ^ {^4}, the inequality (fT3jl 
is a direct consequence of the identity (fT^ . We note that the Gaussian distribution is 
not essential for the proof of the above inequality (fT3j) . 
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The proof of the second inequahty (jHI) is similarly carried out. Derivative by the 
parameter is expressed as in equation by 



-{{Sc)] 



1 d 



{{Sc)] 



dXB (Tb OPb 

Derivative by /^^ is easily calculated as 



d 



NL 



BO 



-{{Sc)] 



NL 



d 



dp. 



B 



i{Sc)] = {JB{SBSc)]-{JB{SB){Sc)]- 



Substitution of equation (fTUj) into equation yields 



1 d 



(Tb dp. 



B 



■{{Sc)] 



J, 



BO 



ctb 



{{SbSc)-{Sb){Sc)] 



- 2aBPB{{SB){SBSc) - {Sb)^{Sc)] 



and the identity (flTjl gives 
d 



ctb- 



' dJBO 

From equation (fT^ we obtain 
d 



{{Sc)] = (rBPB{{SBSc) - {Sb){Sc)] • 



(19) 



(20) 



(21) 



(22) 



da; 



[{Sc)] = 2xb{{SbSc) - {Sb){Sc) - {Sb){SbSc) + {Sb)'{Sc)] • (23) 



B 



The gauge theory yields the following identities on the NL under the same boundary 
conditions as the identity (fTH|) [HI E] : 

{{SbSc)] = {{SbSc)^] 

{{Sb){Sc)] = {{Sb){SbSc)] = {{Sb){Sc){SbSc)] (24) 
{{Sb)'{Sc)] = {{Sb)'{Sc)'] . 
Therefore, we obtain from equation 
d 



dXB 



[{Sc)] = 2xb{{{SbSc) - {SB){Sc)r] > . 



In general, the following relation is proved similarly: 



-{{Sc)'' 



2fc-li 



-{{Sc)''' 



(25) 



(26) 



da;^ da;^ 
= 2A:(2A; - 1)xb{{Sb)''-'{{SbSc) - {Sb){Sc)Y] > . 

It is also possible to prove concavity of the pressure from the derivative of equation 
since the second derivative of the pressure is calculated as 

d 



Xb 



dxc 



{{Sb)] > 



dxBdxc 

where we use equations (jH)) and (fT3j). 



[{Sb)] + 1 + xb-—[{Sb)] > 

dXB 



{B = C) 



(27) 



Letter to the Editor 



5 



3. Discussions on physical consequences 

The inequahty (fTH|l implies that a correlation function [(5^)] is non-negative and 
increases toward unity when the parameter xb of the corresponding subset B tends 
to infinity. This result is reasonable: for large xb, the interaction Jb is almost 
ferromagnetic and the local temperature is nearly zero, and therefore all the spins in 
the subset B become parallel to each other. 

An immediate result from the second inequality is that an arbitrary ra-point 
correlation function, including order parameters, is an increasing function of a parameter 
for any subset. In particular, a two-point correlation function increases with x (which 
is proportional to the inverse temperature and the centre of distribution) and thus the 
correlation length becomes larger as x increases, a natural result. 

The two inequalities have profound consequences on the structure and existence of 
the thermodynamic limit Q y C for both pressure and correlation functions. The first 
inequality ((Zj) tells us that the pressure is monotonically increasing with any xa- Since 
on the NL the Boltzmann weight admits the representation 

eEAsn^AjASA ^ (28) 

where Ja = Ja/o'a is a Gaussian variable with mean [Ja] = xa ^ and variance 
[(jA — xa)'^] = 1, the pressure can be expressed as 

P= [log^e^^e""-^^-^^-* ] , (29) 
s 

and the correlation as 



[{Sc)] 



(30) 



where both P and [{Sc)] are functions only of x's. From equation we see that the 
Xa tune the interaction of the set of spins A. 

We can show the existence of the thermodynamic limit for the pressure per spin 
under free boundary condition. Take a large cube Q (consider for simplicity d = 2) and 
cut it in four identical cubes Qi, i = 1,...,4. We can tune Oj^the interactions among 
the cubes simply by setting all the x's among them equal to zero. In this case the total 
pressure is Pq^ + Pq^ + Pq^ + Pq^. Then we can tune on the interactions from zero to 
their original value and the total pressure becomes Pq. The monotonicity property (|7jl 
gives 

Pn, + Pn, + Pn, + Pn, < Pn , (31) 

which is the well known sub-additivity property for the pressure and implies the existence 
of its density in the thermodynamic limit under two natural assumptions of invariance 
by translation of the Gaussian distributions and the stability boundedness [5]. 

Moreover the second inequality (jH)) can be used to prove the existence of the 
thermodynamic limit for the correlation functions in the same way as the first is used to 
prove the existence of the pressure. In fact, from equation pHjl . we see that when xa = 
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there is no interaction among the subset A, and to increase xa from zero is equivalent 
to add an interaction to the subset A. From the second inequahty (jSJ, addition of an 
interaction for any subset increases all correlation functions. 

This result can be used for the proof that correlation functions have a 
thermodynamic limit under free boundary conditions. Let us consider two finite sets 
Q' C Q. The subset fl is obtained from Q' by adding interactions. Then it is clear from 
the previous argument that 

[{SB)tr^ < [{SB)]t'\ (32) 
This implies that the correlation functions monotonically increase with the system size. 
Since correlation functions are bounded by unity, each of them tends to its unique limit 
as the volume of fl tends to infinity. Therefore the thermodynamic limit for correlation 
functions exists on the NL. 

Monotonicity of the correlation functions with system size can also be proved under 
fixed boundary condition. Fixed boundary conditions {Si = +1 for all boundary sites 
i G dfl) is represented by applying very strong positive magnetic fields to all sites outside 
Q. This is equivalent to xa —>■ oo for A = {i}, i G dQ. Since the set f2 D fi' is obtained 
from Q' by reducing magnetic fields for sites i & we have 

[{Sb)]^^ > [{SB)]t^ (33) 
from the inequality (jHl). Because the inequality (fT^ yields lower bounds for [{Sb}], there 
should exist a well-defined limit as Vl y C The thermodynamic limit of [{Sb)]^^^^ may 
not necessarily be equal to that of [{Sb)]^^'^^'^\ but the former cannot be less than the 
latter. 

The second inequality can also be used for the argument about the location of 
the multicritical points, which are believed to lie on the NL jHl EI; for various lattices. 
For example, let us consider three two-dimensional lattices, the triangular (TR), square 
(SQ) and hexagonal (HEX) lattices. The triangular lattice is obtained from the square 
lattice with addition of bonds, and the square lattice is obtained from the hexagonal. 
Thus, the magnetizations of three lattices satisfy the following relation: 

[(5.)]hex < [(5.)]sQ < . (34) 

Since the multicritical temperature is defined as 

T, = sup{T; [{S,)] > 0} , (35) 

we obtain 

yHEX < ySQ < yTR _ (3g) 

Similarly, the multicritical temperature of the simple cubic lattice is higher then that 
of the square lattice because the former lattice is constructed from the latter by adding 
interactions. 

It is an interesting future problem to prove our results for other models, for example, 
the ± J Ising model, and to develop similar analyses away from the NL. 

After submission of the manuscript we learnt that Kitatani [H] had discussed a 
related problem. 
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Appendix 



In this appendix, we derive equations (fTUjl and (|TT|l using some properties of the Gaussian 
distribution 



Pa{Ja 



exp 



{Ja — JaoY 
2a\ 



Let us consider the following quantity, 
d 



-0{{Ja}) 



B 



I n <^JAPA{JA)jr^O{{JA}) 



(A.i; 



(A.2) 



where 0{{Ja}) is a function of interactions {Ja}- Since Gaussian distribution Pa{,Ja) 
decays rapidly as \ Ja\ — >■ oo^ we may rewrite the right-hand side using integration by 
parts as 



d 



dJ, 



B 



■0{{Ja}) 



I n <^JAPA{JA)dJBPB{JB)TprO{{JA}) (A.3) 



Acq. 
A^B 



I X{^JAPA{JA)dJB{-^PB{JB^O{{JA}) . (A.4) 



At^B 



Calculating the derivative of the Gaussian yields 



d 



Pb{J. 



Jbo — J. 



B 



This relation shows that 
d 



B 



PbUb] 



(9J, 



B 



0{{Ja}) 



n dJAPAiJAy^^^O{{JA}) 



Acn 



a 



B 



-^{[JbOUJa})]- JboPUJa})]} 



(A.5) 

(A.6) 
(A.7) 



which is equation (fTUj). 

For the proof of (fTTj) . we use the property that the Gaussian distribution Pb{Jb) 
is a function of Jb — Jbo, 



_d_ 



B 



■PbU 



B 



d 
dJ 



-Pb{Jb) 



BO 



(A.8) 
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Substitution of (TOIl into yields 



d 



(9J, 



B 



■0{{Ja}) 



[ n dJAPAiJA)dJB (^PBiJB)]0i{JA}) 



Acn 

A^B 

d 



dJBO 

d 



dJ. 



BO 



n dJAPA{JA)0{{JA}) 

Acn 



{OUJa})] . 



(A.9) 
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